Macroscopic dynamics of a Bose-Einstein condensate containing a vortex lattice 
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Starting from the equations of rotational hydrodynamics we study the macroscopic behaviour of a 
trapped Bose-Einstein condensate containing a large number of vortices. The stationary configura- 
tions of the system, the frequencies of the collective excitations and the expansion of the condensate 
are investigated as a function of the angular velocity of the vortex lattice. The time evolution of 
the condensate and of the lattice geometry induced by a sudden deformation of the trap is also 
discussed and compared with the recent experimental results of P. Engels el al, Phys. Rev. Lett. 
89, 100403 (2002). 
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After their experimental realization Jl], 0J quantized 
vortices have become a popular subjet of research in the 
field of ultra cold gases (see for example ||). More re- 
cently, special emphasis has been given to the study of 
configurations containing a large number of vortical lines 
(vortex lattices) j|, || . An important motivation is due to 
the possibility of reaching critical regimes of high angu- 
lar velocities where the degeneracy of the single particle 
levels gives rise to new quantum configurations, analog 
to the quantum Hall effect §, @, |, §. 

The purpose of the present work is to describe the 
macroscopic behaviour of a Bose-Einstein condensate 
containing a vortex lattice in the so called Thomas- Fermi 
regime where the interaction energy is much larger 
than the trapping oscillator energies. A complementary 
investigation concerns the dynamic behaviour of vortex 
lines, where interactions play a role at a more microscopic 
scale @ [|. 

In the presence of a large number of vortex lines a 
useful concept is the so called diffused vorticity given by 



V A v = 2 $7 



(1) 



where v is the velocity field of the fluid and £1 is assumed 
to be uniform. In terms of the density n v of the vortex 
lines per unit surface one has n v — 2QM/h, where M 
is the atomic mass. The quantity $7 corresponds to the 
angular velocity of the sample. The concept of diffused 
vorticity is adequate to describe the dynamics at macro- 
scopic distances, larger than the average distance n v 
between vortices. 

The macroscopic description is provided by the equa- 
tions of rotational hydrodynamics 
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written in the laboratory frame, where n(r, t) is the spa- 
tial density, v(r, t) is the velocity field, V ex t is the external 



potential and g = Ai:h 2 a/M is the coupling costant ex- 
pressed in terms of the s-wave scattering length a. These 
equations generalize the usual equations of irrotational 
hydrodynamics (V A v = 0), which have been success- 
fully applied to study the collective oscillations of the 
condensates in the absence of vortex lines JlJ] . It is worth 
noticing that the hydrodynamic theory is well suited to 
study also non linear effects, including the expansion of 
the gas after release of the confining trap. 

For an axi-symmetric harmonic potential V cxt ~ 
M [uj 2 (x 2 + y 2 ) + lo 2 z 2 }/2 the equilibrium solutions of 
Eqs. (0) and (0) have the form vo = Oq A r and 



n (r) =-•{// 
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where the effective frequency given by ui^ = lo\ — Oq 
originates from the centrifugal effect and 



M = Mo 
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is the chemical potential. In the previous equation 
/jo = (^Who/2) (15Na/aho) 2 ^ 5 is the chemical potential 
in the absence of rotation, uj^o = {i^^i^z) 1 ^ 3 is the av- 
erage oscillator frequency, aho = ^/ft/Afwho is the cor- 
responding oscillator length and N is number of atoms. 
These stationary solutions are defined only for ilo < ■ 
The density profile (|J) exhibits a typical bulge effect pro- 
duced by the centrifugal force, giving rise to the angular 
velocity dependence 



R± 



n 2 



(0) 



for the aspect ratio, where R± and R z are, respectively, 
the Thomas-Fermi radii of the atomic cloud in the radial 
and axial directions. Result (||) can be used to deduce 
the value of f2o from the in situ measurement of the as- 
pect ratio. In Ref. || values up to £1 /luj_ = 0.95 have 
been obtained. In that experiment the condensate takes 
a pancake form even if the confining trap has a cigar 
shape (to z < uj±). 
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The angular momentum per particle carried by the sys- 
tem is given by 



(l z ) = Am (x 2 + y 2 ) 



n 
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and becomes larger and larger as fio — > wx- 

The collective oscillations of the condensate are ob- 
tained by looking for the linearized solutions of Eqs. @) 
and (||) in the form n = uq + Sn and v = vo + Sv with 
Sn and <5v ~ e~ luJt . 

If one chooses 5n — a±(x ± iy) 2 , <$v = ay-V(;r ± iy) 2 , 
corresponding to quadrupole excitations with m = ±2 
where m is the third component of angular momentum, 
one easily finds the dispersion law |2ll 



u(m = ±2) = \ 2uo]_ - fi 2 ± fi 



(8) 



These frequencies have been directly measured in the ex- 
periment of confirming with high accuracy the pre- 
dictions of theory. The splitting 2 fio between the two 
frequencies agrees with the sum rule result |lq] 



u){m = +2) — Lu(m — —2) = — — 
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as can be seen using the rigid body expression (^) for the 
angular momentum. 

When fio — > u± one finds u)(m = +2) — > 2ljj_, while 
u)(m = —2) — > 0, reflecting the tendency of the system to 
become unstable against static quadrupole deformations. 

In addition to the m = ±2 oscillations, it is interesting 
to discuss the behaviour of the m = modes [^2) . These 
result from the coupling between the axial and radial 
motion of the condensate. By looking for solutions of 
the form 5n(r) = ciq + a±(x 2 + y 2 ) + a z z 2 and 5v(r) = 
Stt A r + V[aj_(x 2 + y 2 ) + a z z 2 ] one finds that the two 
decoupled frequencies are given by 
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When fio = one recovers the solutions of Q|, while 
for £Iq — > lo±_ the two frequencies approach the value 
2uj± (radial compressional mode) and y/3oj z (axial com- 
pressional mode). The latter value coincides with the 
frequency predicted by theory for a pancake geometry 
(u) z ^> u±) in the absence of vortices. 

The tendency of the system to become unstable against 
a quadrupole deformation when fio — > lo± has been used 
in the recent experiment of |l5[ ] to induce large deforma- 
tions in the condensate and to explore the consequences 
on the geometry of the vortex lattice. 

In the presence of a static deformation of the trap 
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FIG. 1: Deformation S — a/fig of the cloud for the steady- 
state solutions of Eqs. (|^) and (g|) as a function of the angular 
velocity flo of the vortex lattice, for a fixed value e = 0.036 
of the trap deformation. The dotted line is the asymptotic 
solution for e — > 0. 



the stationary solutions of the hydrodynamic equations 
(||) and (||) in the laboratory frame change in a quite 
interesting way. On the one hand the velocity field takes 
the form 



v (r) = n Ar + aV(xy) , 



(12) 



containing a crucial irrotational component fixed by the 
parameter a. On the other hand the equilibrium profile 
assumes the deformed shape 
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with the effective frequencies defined by w 2 „, = u 2 



a 2 — and the chemical potential given by [i = 
/j,o(cu x LOy/uj x ujy) 2 / 5 . Use of the equation of continuity 
yields the following third order equation for a 



+ ol{uj\ - - euj']_ fi = 



(14) 



whose solution, expressed in terms of the deformation of 
the condensate 



(y 2 - x 2 ) 

(x 2 + y 2 ) 



a 
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is reported in Fig. [TJ for the value e = 0.036. The other 
solutions of the third order equation (|l4|) should be ex- 
cluded because they correspond to negative values of 



~j 2 



only irrotational flow was considered and more stationary 
solutions were found to be available in the rotating frame. 
For e ^ 0, Fig. |l| shows that, differently from the axisym- 
metric case, a deformed system can support values of fio 
larger than u;j_j the irrotational term of Eq. ( |l2| ) pro- 
viding the crucial compensation to the centrifugal effect 
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generated by the rotational component. When e — > and 



> w ± one finds u£ 



and the branch approaches 



the asymptotic solution S = \J\ — (a;j_/^o) 2 - 

In the experiment of pl| the static deformation (|l] 
was switched on suddenly at some initial time t — 0. This 
produces a time tependent perturbation which drives the 
system far from the initial axi-symmetric configuration. 
In the following we will discuss the prediction of Eqs. @j 
and (H) for the time evolution of the condensate shape 
as well as of the vortex lines which are assumed to follow 
the motion of the fluid. The equations are easily solved 
in the non linear regime by looking for the ansatz 



(16) 



n = cio + a x x 2 + a y y 2 + a z z 2 + a xy xy 
v = fi A r + "V(a x x 2 + a y y 2 + a z z 2 + a xy xy), (17) 

where now f2, as well as the coefficients Oj and ctj, is time 
dependent. In Fig. ^(a) we report the predictions for the 
angle between the principal axis of the condensate and 
the y axis of the trap, with the initial value of the angular 
velocity Slo = 0.95 ui±. The angle first increases fast fol- 
lowing the direction of the vortex flow. This happens be- 
cause for short times the perturbation excites with equal 
strength both the rn = +2 and m = —2 modes and 
hence the system exhibits a precession with angular ve- 
locity proportional to Lu(m = +2) — uj(m = —2). For 
longer times the response is governed by the static ef- 
fect and is dominated by the m = — 2 mode, producing 
a precession in the opposite direction whose period be- 
comes longer and longer as the value of e decreases. The 
deformation of the condensate as a function of time is 
plotted in Fig. |^(b) for the same conditions: rather large 
values are reached, in agreement with the experimental 
findings of ]l5|. For example at t ~ 2.5 T± the authors 
of E5l report the value 5 ~ 0.4. The coupling with the 
fast m = +2 mode is evident in both figures and results 
in the small oscillations of higher frequency. 

Once the solutions of the hydrodynamic equations are 
known one can also study the time evolution of the vor- 
tex lines wich are assumed to move following the velocity 
field v (Fig. ||). We have initially assumed [Fig. ||(a)] an 
hexagonal geometry (Abrikosov lattice |l9|). In the ab- 
sence of deformation the vortex lines rotate rigidly with 
angular velocity Qq. However, once the deformation of 
the condensate takes place, the vortex lattice changes its 
geometry since it is forced to follow the stream lines of 
the velocity flow (17) wich contains an irrotational com- 
ponent. For example, in Fig. ||(b) one sees the formation 
of a near orthorhombic lattice, while for large deforma- 
tions one observes peculiar structures characterized by 
the compression of vortex lines along the short radius 
[Fig. ||(c)]. In our description the vortex lines are just 
points of the fluid moving with velocity v, so we do not 
account for the possible deformation of the density pro- 
files induced by the vicinity of the vortex lines. These 
effects are expected to give rise to a significant suppres- 
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FIG. 2: Time dependence of the condensate after switching 
on the trap deformation with e — 0.036: (a) evolution of the 
angle 9(t) between the principal axis of the cloud and of the 
trap in the x — y plane; (b) evolution of the cloud deformation 
8(t). The axial frequency of the trap is uj z — 0.65 u>±, while 
the initial configuration is the stationary solution for e = 
and f2 = 0.95 uj±. Time is measured in units of Tj_ = 2n/u}±. 
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FIG. 3: Evolution of the vortex lattice after switching on 
the trap deformation e = 0.036: (a) t = 0, (b) t = 2.8 T ± , 
(c) t — 3.9 T±. The peculiar geometry shown in (c) appears 
when the cloud deformation 5 is large and one of the lattice 
vector is aligned along the short axis of the condensate. The 
trap parameters are the same as in Fig. [| 



sion of the density in the region between two vortices 
when their relative distance becomes small, with the ap- 
pearence of stripe like configurations as observed in the 
experiments of []l5| and recently discussed in [|l3) . 

An important role in the experimental detection of vor- 
tices is played by the expansion of the condensate. The 
vortex cores are indeed too small to be imaged in situ. 
It is therefore important to look at the behaviour of the 
condensate after the sudden release of the trap. This 
can be investigated by solving the hydrodynamic equa- 
tions Q) and (||) setting V cxt — at the release time. 
The general solution is still given by the form (|l^) and 
( jl7| ) . We have first studied the time dependence of the 
aspect ratio (^), starting from an axisymmetric configu- 
ration (e = 0). Using the trapping parameters of Rcf. pl| 
we find that R±/R z increases in time (see Fig. |j). This 
dramatically differs from the expansion of a non rotat- 
ing condensate, which would instead transform a pan- 
cake cloud into a cigar shape one. This behaviour is the 
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FIG. 4: Time evolution of the aspect ratio (^[) after re- 
leasing the axisymmetric trap for a condensate rotating at 
fio = 0.95 uij_ (solid line) and Qo ~ 0.90 u)± (dashed line). 
The anisotropy of the trap is A = uj z /lj± = 0.65. 



consequence of the large radial rotational kinetic energy, 
which produces a fast expansion in the radial direction, 
and becomes more and more pronounced as Qq ~~ ► 
We have also investigated the time evolution of the cloud 
deformation S(t) in the x — y plane, by switching off the 
confining potential when the perturbation ( |TT|) has pro- 
duced a sizable condensate deformation [see Fig. §(b)]. 
Using the experimetal parameters of flfi ] we find that 
the value of the condensate deformation remains almost 
constant during the expansion. This is again the conse- 
quence of the centrifugal effect, which is stronger in the 
direction of the long axis and compensates the effects of 
the pressure gradient. 

Let us finally discuss the conditions of applicability 
of the Thomas-Fermi approximation used in the present 
work. If £}q becomes too close to u>± the condition 
fi fiw± , ftto z required to apply the Thomas- Fermi ap- 
proximation is no longer valid [see Eq. (^)]. In partic- 
ular if [i ~ huj± the size of the vortex core, fixed by 
the healing length £ = %JsJ2M\i, becomes comparable 
to the average distance between vortices, fixed by the so 
called effective magnetic length £q = y/H/2MClo. In 
fact one has (,/l& — \/MIq/i1 ~ \J kuj±^ / \i. In the exper- 
imental conditions of Rcf. Jig] one has /io ~ 27 fowho and 
the Thomas-Fermi approximation is well satisfied also 
for f^o = 0.95 Deviations from the Thomas-Fermi 
regime require smaller values of \i§ and/or larger values 
of f^o / <^j_ and should show up in a different behaviour of 
the collecitve frequencies. 

Let us also mention that in order to observe quantum 
Hall features one should reach even more extreme con- 
ditions, where the filling factor v = N/N v , given by the 
ratio between the number of atoms and the number of 
vortices, is of the order of unity pSj]. Under these condi- 



tions the low frequency oscillations in the rotating frame 
are expected to exhibit a chiral behaviour p0[ , similar 
to what happens to the edge excitations in the quantum 
Hall effect. 

In conclusion, we have studied the stationary configu- 
rations and the dynamics of a Bose-Einstein condensate 
both in an axisymmetric and in an anisotropic static trap, 
using the equations of rotational hydrodynamic and the 
concept of diffused vorticity. Special emphasis has been 
given to the collective oscillations, wich turn out to be 
significantly affected by the rotation of the cloud. In ad- 
dition, our macroscopic description gives valuable insight 
on the evolution of the vortex lines recently observed in 
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Eq. (^) can be generalized to excitations of higher 
multipolarity £ and m = ±£. One finds jTj| lu± = 
■s/£u}\ - {£-1)0,1 ± V ~ !) Q o- 

The m = ±1 modes, wich involve a dynamic coupling 
with the rotation of the axis of vorticity (tilting mode) , 
will be discussed elsewhere |L7| . 

For example, assuming a 3-D Thomas-Fermi profile, one 
finds v = (2/15)(p/hQo)(Rz/a) and hence v 3> 1 even if 
H ~ hO,a. 



